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Abs tract 

The results of an experimental and analytical 
study of the effects of initial imperfections on the buck- 
ling and postbuckling response of unstiffened thin- 
walled compression-loaded graphite-epoxy cylindrical 
shells are presented. The shells considered in the study 
have four different shell- wall laminates and two differ- 
ent shell-radius-to-thickness ratios. The shell-wall lam- 
inates include two different orthotropic laminates and 
two different quasi-isotropic laminates. The shell-radi- 
us-to-thickness ratios include shell-radius-to-thickness 
ratios equal to 100 and 200. The results identify the ef- 
fects of traditional and nontraditional initial imperfec- 
tions on the nonlinear response characteristics and 
buckling loads of the shells. The traditional imperfec- 
tions include the geometric shell-wall mid-surface im- 
perfections that are commonly discussed in the 
literature on thin shell buckling. The nontraditional im- 
perfections include shell-wall thickness variations, lo- 
cal shell-wall ply-gaps associated with the fabrication 
process, shell-end geometric imperfections, nonuni- 
form applied end loads, and variations in the boundary 
conditions including the effects of elastic boundary 
conditions. A high-fidelity nonlinear shell analysis pro- 
cedure that accurately accounts for the effects of these 
traditional and nontraditional imperfections on the non- 
linear response characteristics and buckling loads of the 
shells is described. The analysis procedure includes a 
nonlinear static analysis that predicts the stable re- 
sponse characteristics of the shells, and a nonlinear 
transient analysis that predicts the unstable response 
characteristics. The results of a local shell-wall stress 
analysis used to estimate failure stresses are also de- 
scribed. 

Ij uro d uc tioa 

The increasing need to produce lighter-weight 
aerospace shell structures has led to the use of advanced 
material systems in new structural designs, and im- 
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proved design methods appropriate for these advanced 
material systems are needed. The high strength-to- 
weight and high stiffness-to- weight ratios of advanced 
composite materials offer significant weight reduction 
potential for aerospace structures. In addition, the use 
of advanced composite materials allows the designer to 
tailor the stiffness properties of composite structures to 
obtain structurally efficient designs. Designers often 
use a design-level analysis procedure with empirical 
data to develop new structural designs for strength and 
buckling-critical structures. The traditional approach 
for designing thin-walled buckling-resistant isotropic 
shell structures is to predict the buckling load of the 
shell with a deterministic analysis, and then to reduce 
this predicted load with an empirical “knockdown’' fac- 
tor (e.g., Ref. 1). The empirical knockdown factor is in- 
tended to account for the difference between the 
predicted buckling load and the actual buckling load for 
the shell determined from tests. A linear bifurcation 
buckling analysis is often used for the design-level 
analysis, and this analysis is usually based on nominal 
structural dimensions and material properties of an ide- 
alized, geometrically perfect shell. The design knock- 
down factor used in the design of buckling-resistant 
shells is often based on the “lower bound” design rec- 
ommendations reported in Ref. 1 . This design philoso- 
phy can result in overly conservative designs for these 
structures, and it can potentially even result in uncon- 
servative designs if the empirical data are not represen- 
tative of the design of interest. While it is generally 
recognized that initial geometric shell-wall imperfec- 
tions are a major contributor to the discrepancy between 
the predicted shell buckling loads and the experimental- 
ly measured shell buckling loads (e.g., Refs. 2-6), the 
traditional sources of design knockdown factors do not 
include data or information related to the sensitivity of 
the response of a shell to various forms of imperfec- 
tions. In addition, the traditional sources of design 
knockdown factors for predicting shell buckling loads 
do not include information for shell structures made 
from advanced composite materials. Recent studies 
(e.g., Refs. 7-11) have shown that traditional initial geo- 
metric shell-wall imperfections, and other nontradition- 
al forms of imperfections or variations in geometric and 
material parameters, loading conditions, and boundary 
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conditions can significantly affect the buckling load of 
a compression-loaded composite shell structure. Other 
recent studies (e.g., Ref. 12) have identified the effects 
of fabrication related anomalies on the buckling loads 
of metallic shells. The effects of these traditional and 
nontraditional classes of initial imperfections on the 
buckling of composite shells are generally not well un- 
derstood by structural engineers and designers. Modern 
high-fidelity nonlinear analysis procedures (e.g.. Ref. 
13) offer the opportunity to improve some of the engi- 
neering approximations that are used in the design and 
analysis of shell structures, and to provide insight into 
the effects of traditional and nontraditional imperfec- 
tions on the response of compression-loaded composite 
shell structures. 

The present paper describes the results of an ex- 
perimental and analytical study of the effects of tradi- 
tional initial geometric shell-wall imperfections, and 
the effects of nontraditional initial imperfections and 
variations in other nontraditional geometric and materi- 
al parameters, loading conditions, and boundary condi- 
tions on the buckling response of unstiffened thin- 
walled compression-loaded graphite-epoxy cylindrical 
shells. The results of four graphite-epoxy shells with 
different shell-wall laminates are presented. The shell- 
wall laminates include two different orthotropic lami- 
nates and two quasi-isotropic laminates. The two 
orthotropic shells and one of the quasi-isotropic shells 
have shell-radius-to-thickness ratios equal to 200, and 
the other quasi-isotropic shell has a shell-radius-to- 
thickness ratio equal to 100. Traditional shell- wall geo- 
metric imperfections and several nontraditional imper- 
fections were measured, and representations of these 
imperfections have been included in nonlinear analyses 
of the shells. These nonlinear analyses were conducted 
with the geometrically nonlinear STAGS finite element 
analysis code (Ref. 13). The effects of initial geometric 
shell-wall imperfections, shell-wall thickness varia- 
tions, shell-end geometric imperfections, nonuniform 
applied end loads, and variations in the boundary con- 
ditions, including the effects of elastic boundary condi- 
tions, on the buckling response of these thin-walled 
composite shells are discussed in the present paper. 
The effects of fabrication anomalies caused by small lo- 
cal gaps between adjacent pieces of the graphite-epoxy 
material in some of the shell-wall laminate plies are 
also discussed. The response of a quasi-isotropic shell 
with a shell-radius-to-thickness ratio equal to 100 is 
presented as an example of a thin-walled shell that has 
material failures before it buckles, and a local stress 
analysis of the shell wall is used to show that these local 
material failures can result in the overall failure of the 
shell. The results of the study are used to illustrate the 
significance of these nontraditional initial imperfec- 


tions on composite shell response characteristics. The 
nonlinear shell analysis procedure used to predict the 
nonlinear response and buckling loads of the shells is 
described, and the analysis results are compared with 
the experimental results. The use of this nonlinear shell 
analysis procedure for determining accurate, high-fi- 
delity design knockdown factors for shell buckling and 
collapse, and for determining the effects of variations 
and uncertainties in shell geometric and material pa- 
rameters on shell buckling loads is discussed. 

Test Specimens, Imperfection Measurements, and Test 
Apparatus and Tests 

Test Specimens 

The specimens tested in this investigation were 
fabricated from 0.005-in. -thick AS4/3502 graphite-ep- 
oxy preimpregnated unidirectional tape material made 
by Hercules, Inc. The nominal unidirectional lamina 
properties of a typical 0.005-in. -thick ply with a fiber 
volume fraction of 0.62 are as follows: longitudinal 
modulus Ei = 19.5 Msi, transverse modulus E 2 = 1.45 
Msi, in-plane shear modulus Gy 2 = 0.813 Msi, and ma- 
jor Poisson’s ratio v 12 = 0.30. The material was laid up 
on a 15.75-in. -diameter mandrel and cured in an auto- 
clave to form four shells with different shell-wall lami- 
nates including an 8-ply axially stiff [+45/0 2 ] s 
laminate, an 8-ply circumferentially stiff [+45/90 2 ] s 
laminate, an 8-ply quasi-isotropic [+45/0/90] s lami- 
nate, and a 16-ply quasi-isotropic [+45/0/90] 2s lami- 
nate. The resulting four shells are referred to herein as 
shells or specimens Cl, C2, C3 and C4, respectively. 
These specimens had a nominal length of 16.0 in. and a 
nominal radius of 8.0 in. The 8-ply and 16-ply speci- 
mens had a nominal shell-wall thickness of 0.04 in. and 
0.08 in., respectively, and had shell-radius-to-thickness 
ratios of 200 and 100, respectively. Both ends of the 
specimens were potted in an aluminum-filled epoxy 
resin to assure that the ends of the specimens did not fail 
prematurely during the test. The potting material ex- 
tended approximately 1.0 inch along the length of the 
specimens at each end resulting in a test section that 
was approximately 14.0 in. long. The ends of the spec- 
imens were machined flat and parallel to assure proper 
load introduction during the tests. A photograph of a 
typical specimen and the specimen coordinate system 
used to represent the corresponding geometry is shown 
in Fig.l. The shell length, test-section length, radius, 
and thickness are designated as L, Lj, R and r, respec- 
tively. 

Imperfection Measurements 

Three-dimensional surveys of the inner and outer 
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shell-wall surfaces of the specimens were made prior to 
testing the specimens to determine their initial geomet- 
ric shell-wall imperfection shapes and shell-wall thick- 
ness distributions. Measurements were taken over a 
uniform grid with increments of 0.125 in. in the axial 

direction and 0.139 in. (approximately 1° of arc) in the 
circumferential direction over the exposed surfaces of 
the specimens. The inner surface measurement was 
used to determine the initial geometric shell-wall im- 
perfection shape of a specimen, and the difference be- 
tween the outer and inner surface measurements was 
used to determine the shell-wall thickness distribution, 
A contour plot of the nondimensionalized initial geo- 
metric shell-wall mid-surface imperfections w o (;c,0) 

for specimen C3 is shown in Fig. 2. The measured 
shell-wall imperfection w 0 is nondimensionalized by 
the average measured shell-wall thickness t ave = 
0.0381 inches. These results indicate that the initial 
geometric shell-wall imperfection is periodic in the cir- 
cumferential direction and has slight deviations in the 
axial direction. The amplitude of the imperfection var- 
ies from + \34\t ave to -1.535 t ave . A contour plot of 
the nondimensionalized shell-wall thickness variation 
r o (jc,0) for specimen C3 is shown in Fig. 3, where the 
measured thickness value t 0 is nondimensionalized by 
the average measured shell-wall thickness t ave . These 
results indicate that the shell-wall thickness, and hence 
the laminate stiffnesses, varies significantly over a short 
distance. The thickness varies from 0.928 to 1.321 
times t ave . Most of the thickness variation is attributed 
to local variations in the resin content of the laminate 
associated with the fabrication process. However, the 
darker angular pattern in the thickness distribution is at- 
tributed to small gaps between adjacent pieces of graph- 
ite-epoxy tape in some of the laminate plies that were 
generated during the lay-up and curing processes. Such 
a region is referred to herein as a lamina ply-gap or a 
ply-gap. In the particular case of shell C3, such a local- 
ly thin region in the shell wall consists of a 7-ply-thick 
laminate rather than the nominal 8-ply-thick laminate. 
For the case where one ply-gap intersects another ply- 
gap, the shell wall consists of a 6-ply-thick laminate. 
Moreover, these locally thin shell-wall regions have a 
significant shell-wall mid-surface eccentricity, and 
have reduced stiffnesses relative to the rest of the shell 
wall. Typically, a circumferential ply-gap is caused by 
a gap between two adjacent 90° pieces of tape in the 
ply, a 45° or helical ply-gap is caused by a gap between 
two adjacent 45° pieces of tape, and an axial ply-gap is 
caused by a gap between two adjacent 0° pieces of tape. 


Magnified cross-sectional views of typical ply-gaps in 

a -45° outer ply of a specimen laminate and in a 90° 
mid-surface ply are shown in Figs. 4a and 4b, respec- 
tively. The widths of the ply-gaps shown in Figs. 4a 
and 4b are on the order of 0.02 in. or approximately 
equal to half the shell-wall thickness, and the ply-gap 
depth is approximately 0.005 in. or approximately 
equal to the nominal lamina ply thickness. In addition, 
there often exist regions adjacent to the ply-gaps in 
which a gradual thickness reduction of the lamina ply 
occurs as shown in Fig. 4a. These regions of lamina ply 
thickness reduction can extend an additional 0.02 to 
0.04 in. from the edges of the ply-gap. Lamina ply-gaps 
with gap widths as large as 0. 1 in. have been observed 
in some of the shell specimens. The lighter angular pat- 
terns in the thickness contour plot are caused by locally 
thickened regions of the outermost plies of the laminate 
that develop during the curing process to form outer 
shell-wall surface ridges. A magnified cross-sectional 
view of such a region is shown in Fig. 4c. Some of the 
shell-wall thickness features, such as lamina ply-gaps, 
are smaller than the imperfection-measurement grid 
spacing used in this study, and as a result, some of the 
smaller thickness variation features may not have been 
included in the measurements. 

Measurements of the top and bottom loading sur- 
faces of the specimens were made every degree around 
the circumference of the specimens to determine the 
variation in the shell-end or loading-surface geometry. 
Typical top and bottom shell-end or loading-surface ge- 
ometry variations for specimen C3 are denoted by 
S top (Q) and 5^(0), respectively, and are shown in Fig. 
5. The maximum amplitude of this shell-end or load- 
ing-surface variation is approximately 0.0015 inches, 
which is approximately 4% of t ave or 0.01% of the 
specimen length. 

Test Apparatus and Tests 

The specimens were instrumented with electrical 
resistance strain gages, and direct-current differential 
transducers (DCDT’s) were used to measure displace- 
ments. Three non-collinear DCDT’s were positioned at 
three corners of the upper loading platen of the test ma- 
chine and used to measure the end-shortening displace- 
ment A and the rotations § y and <|>, of the loading 
platen as illustrated in Fig. I . A typical set of measured 
loading platen rotations is shown in Fig. 6. These re- 
sults indicate that significant upper platen rotations oc- 
cur from the onset of loading up to a load value of 
approximately 6,000 lbs. These rotations are attributed 
to an initial misalignment of the upper loading platen 
and the specimen. The rotations of the movable upper 
loading platen reach a steady state at a load value of 
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6,000 lbs., and the loading of the specimen, for the most 
part, continues without additional upper loading paten 
rotations from 6,000 lbs. up to the buckling or collapse 
load. During the collapse event, the upper loading plat- 
en undergoes an additional amount of rotation. The 
shadow moir£ interferometry technique was used to ob- 
serve the shell-wall prebuckling, buckling and post- 
buckling normal (perpendicular to the shell outer 
surface) deformation patterns. All data were recorded 
with a data acquisition system, and the moire patterns 
were recorded photographically and on videotape. 

The specimens were loaded in compression with a 
300,000-lb hydraulic universal-testing machine by ap- 
plying an end-shortening displacement to the shell ends 
or loading surfaces of the specimens. To control the 
load introduction into the specimens, the upper loading 
platen was aligned with the loading surface of the spec- 
imen as well as possible before the test by adjusting lev- 
eling bolts in the corners of the upper loading platen 
until strains measured by selected strain gages on the 
specimens indicated a uniform axial strain distribution 
in the shell wall. The specimens were loaded until gen- 
eral instability or failure of the shells occurred. 

Finite-Element Models and Analyse s 
Nonlinear Analysis Procedure 

The shells considered in this study were analyzed 
with the STAGS (STructural Analysis of General 
Shells) nonlinear shell analysis code. 13 STAGS is a fi- 
nite-element code developed for the nonlinear static 
and dynamic analysis of general shells, and includes the 
effects of geometric and material nonlinearities in the 
analysis. The code uses both the modified and full 
Newton methods for its nonlinear solution algorithms, 
and accounts for large rotations in a shell by using a co- 
rotational algorithm at the element level. The Riks 

pseudo arc-length path-following method 14 is used to 
continue a solution past the limit points of a nonlinear 
response. With this strategy, the incrementally applied 
loading parameter is replaced by an arc-length along the 
solution path, which is then used as the independent 
loading parameter. The arc-length increments are auto- 
matically adjusted by the program as a function of the 
solution behavior. The transient analysis option in 
STAGS uses proportional structural damping and an 
implicit numerical time-integration method developed 

by Park. 15 Additional information on the transient 
analysis procedure can be found in Ref. 16. 

The prebuckling, buckling and postbuckling re- 
sponses of the shells were determined using the follow- 
ing analysis procedure. The prebuckling responses 
were determined using the geometrically nonlinear 
quasi-static analysis capability in STAGS. The Riks 


pseudo arc-length path-following method in STAGS 
was used to compute the initial shell response until just 
before buckling occurred. The unstable buckling re- 
sponse of the shell was predicted using the nonlinear 
transient analysis option of the code. The transient 
analysis was initiated from an unstable equilibrium 
state close to the limit point by incrementing the end 
displacement by a small amount. An initial time step of 
1.0E-8 seconds was used in the analysis, and the time 
step was automatically adjusted by the program as a 
function of the solution behavior. The transient analy- 
sis was continued until the kinetic energy in the shell 
had dissipated to a negligible level, which indicated that 
the transient response had attenuated. Once the tran- 
sient analysis had attenuated to a near-steady-state so- 
lution, the load relaxation option of the code was used 
to establish a static equilibrium state. Conventional lin- 
ear bifurcation buckling analysis results were also de- 
termined with STAGS for comparison with the 
nonlinear response results. 

Finite-Element Models 

A typical finite-element model of a specimen is il- 
lustrated in Fig. I. The standard 410 quadrilateral ele- 
ment from the STAGS element library was used in the 
modeIs._The elements of the finite-element mesh are 
approximately 0.2-in. by 0.2-in. square. Each element 
has four integration points, which are distributed in 
such a way as to provide a modeling resolution of ap- 
proximately 0.1-in. by 0.1 -in square. This integration- 
point spacing is on the order of the measurement- point 
spacing used when measuring the initial geometric im- 
perfections of the specimens. This highly refined mesh 
is necessary to model rapidly varying geometric and 
material parameters such as nonuniform shell-wall 
thicknesses and lamina stiffness properties. A typical 
finite-element model contained approximately 100,000 
degrees of freedom. 

Geometrically perfect and imperfect shells were 
analyzed in the present investigation. Nominal shell 
geometry, laminate thickness, lamina mechanical prop- 
erties, and boundary conditions were used for the finite- 
element models of the geometrically perfect shells. The 
nominal boundary conditions consist of setting the cir- 
cumferential and normal displacements v and w equal 
to zero in the 1.0-in. -long potted boundary regions of 
the shell illustrated in Fig. 1, setting u{LI 2, 0) = 0, and 
applying a uniform end-shortening u{-U 2, 0) = A. 

The geometrically perfect finite-element models 
were modified to include the effects of the measured 
shell imperfections in order to simulate more closely 
the response of the specimens. These modeling modi- 
fications include the effects of the measured initial geo- 
metric shell-wall mid-surface imperfections, shell-wall 
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thickness variations, local shell-wall lamina ply-gaps, 
thickness-adjusted lamina properties, elastic boundary 
conditions, shell-end geometric imperfections, and 
nonuniform end loads. 

The initial geometric shell-wall mid-surface im- 
perfection w„(*,0) is included in the finite-element 
models by introducing an initial normal perturbation to 
each node of the mesh by using a user-written subrou- 
tine with STAGS for that purpose. The user-written 
subroutine uses a linear interpolation algorithm that cal- 
culates the value of the imperfection for the coordinates 
of each finite-element node based on the measured 
shell-wall data. 

The shell-wall thickness r, mid-surface eccentric- 
ity ecZy and lamina material properties E 2 , G,2 
and V /2 are adjusted at each integration point of each 
element in the finite-element models. The shell-wall 
mid-surface eccentricity is calculated relative to the av- 
erage shell-wall mid-surface as illustrated in Fig. 7; that 
is, eczixfi) = -0.5(7^ - t 0 (x y 0)). The lamina properties 
are adjusted by using the rule of mixtures. In the rule- 
of-mixtures calculations, it is assumed that any varia- 
tion in the lamina ply thickness from the nominal thick- 
ness is due to a variation in resin volume only, and that 
the fiber volume remains constant for each ply. How- 
ever, there are several assumptions and approximations 
related to modeling the geometry and stiffnesses of a 
lamina ply-gap detail that were used in the analysis. 
First, the finite-element models are limited to modeling 
the shell-wall thickness variation as discrete step- 
changes, as illustrated in Fig. 7, and the resolution of 
the thickness variation is limited by the finite-element 
integration point spacing (i.e., 0.1 in.). Results from a 
study illustrating the potential effects of these modeling 
approximations and mesh refinement on the response of 
a shell with a ply-gap indicate that this modeling ap- 
proach can affect the predicted response. It was found 
that these modeling approximations can lead to an arti- 
ficial increase in the size of the measured thickness de- 
tails, e.g., the width of a lamina ply-gap, resulting in 
misrepresentation of the local bending stiffnesses and 
the mid-surface eccentricity of the shell wall at that par- 
ticular point. In addition, the mesh refinement and in- 
tegration point spacing used in the present models tend 
to provide models that are overly stiff in bending by a 
small amount when they are used for modeling a ply- 
gap detail. These modeling approximations can have a 
direct influence on the local bending response of the 
shell and can result in as much as a 2% variation in the 
predicted buckling loads. 

The second assumption is related to the modeling 
of the stiffnesses of the lamina ply-gap. Two modeling 
approaches were considered. One approach is to model 


each ply-gap with a local reduction in thickness, as 
measured, including the local eccentricity ecz , and re- 
ducing the appropriate number of lamina plies in the 
shell-wall laminate model, hence, reducing the local 
stiffnesses associated with the local ply-gap detail. The 
second approach is to model the ply-gap with the local 
as-measured reduced shell-wall thickness and the corre- 
sponding local mid-surface eccentricity, but assuming 
that the local thickness reduction is due to a reduction 
in the resin volume only and, consequently, keeping the 
fiber volume constant. This approach neglects the local 
stiffness reduction associated with the ply-gap details. 
The former modeling approach is very time consuming 
to implement since it requires the manual definition of 
each shell-wall laminate associated with each specific 
ply-gap in the finite-element model. The latter model- 
ing approach is much easier to implement in the model 
by using a user-written subroutine compatible with the 
STAGS finite-element analysis code. As a result, the 
latter modeling approach was used in the present study. 
Results from a study of this modeling assumption indi- 
cate that neglecting the local stiffness reduction associ- 
ated with a ply-gap would cause only slight differences 
in the magnitude of the local bending response and no 
more than a 2% variation in the predicted buckling load 
for a compression-loaded shell. In addition, the results 
of this modeling study indicated that the local shell-wall 
mid-surface eccentricity is the most important feature 
of the ply-gap detail for these stability critical prob- 
lems, and this eccentricity was included in the models 
for the results presented herein. 

To provide a better simulation of the constraints 
provided by the potting material at the ends of the spec- 
imens, effective axial and radial potting-support stiff- 
nesses were determined for each shell specimen using a 
two-dimensional generalized plane-strain finite-ele- 
ment analysis of the potting-material-shell-wall detail 
shown in Fig. 8a. Material properties of the potting 
compound were characterized by Iosipescu tests of the 

potting material reported by Weiland et al. 17 The nom- 
inal properties of the potting material are as follows: 
Young’s modulus E= 1.15 Msi, shear modulus Gj 2 - 
0.36 Msi, and Poisson’s ratio v = 0.59. The effective 
laminate axial stiffness E x was used in the shell-wall 
models which include E x equal to 11.83, 3.60, 8.07, 
and 16.14 Msi, for shells Cl, C2, C3 and C4, respec- 
tively. To determine the effective axial and radial pot- 
ting-support stiffnesses, the following numerical 
experiments were conducted. In the first numerical ex- 
periment, a uniform axial displacement A was applied 
to an unpotted end of the shell model, and the resulting 
predicted axial strain response, shown in Fig. 8b, was 
used to calculate an effective axial stiffness K A for the 
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portions of the shell wall supported by the potting ma- 
terial. The results indicate that the axial strain response 
of the shell wall supported by the potting material is a 
function of shell-wall laminate modulus. In particular, 
the decay length of the strain response varies inversely 
with the axial laminate stiffness. The predicted effec- 
tive axial potted-shell stiffnesses are 26.55, 6.73, 26.26 
and 32.37 Msi for shells C 1 , C2, C3 and C4, respective- 
ly. In the second numerical experiment, a unit concen- 
trated force F z was applied to the shell model, and the 
resulting w (radial or normal) displacement and nodal 
force at the junction of the shell wall and the potting 
material was used to calculate an effective radial sup- 
port stiffness K R . The effective radial potting-support 
stiffness was predicted to be approximately equal to 
1.0E5 lbf/in. 

Nonuniform end loading of a specimen is attribut- 
ed to initial specimen-end or loading-surface imperfec- 
tions and to upper loading-platen rotations that are 
measured during the experiment. First, the measured 
upper and lower specimen-end or loading-surface im- 
perfections 5 t0 p(Q) and <5^ f (0), respectively, were in- 
cluded in the finite-element model by introducing an 
initial in-plane axial perturbation to the nodes at the 
loaded ends of the shell. Then, the compression load 
was applied to the shell in two parts. The nonuniform 
specimen-end imperfections, -S top (B) and 
Sb 0[ (&), were applied as displacements to the upper and 
lower ends of the shell, respectively, at the beginning of 
the analysis to simulate a full contact condition between 
the shell ends and the loading platens. Then, the exper- 
imentally measured end-shortening displacement A 
and upper loading-platen rotations § y and <j) z were ap- 
plied to the upper shell end or loading surface while 
holding the lower loading surface fixed as illustrated in 
Fig. 1; that is, w(-L/2, 0) = A + R cos<j) y cos0 + R cos 
<t>z sin 0 -5, op (Q) and u(U2, 0) = -S bol (Q). 

Results and Discussion 

Analytically predicted and experimentally mea- 
sured results for the four compression-loaded graphite- 
epoxy cylindrical shells considered in this study are 
presented in this section. The shell-wall laminates of 
the four shells include two different orthotropic lami- 
nates and two different quasi-isotropic laminates. The 
two 8-ply orthotropic shells, Cl and C2, and the 8-ply 
quasi-isotropic shell, C3, have shell-radius-to-thickness 
ratios equal to 200. The 16-ply quasi-isotropic shell, 
C4, has a shell-radius-to-thickness ratio equal to 100. 
The predicted results were obtained from finite-element 
models of geometrically perfect shells and shells that 
include initial geometric shell-wall mid-surface imper- 
fections, shell-wall thickness variations and thickness- 


adjusted lamina properties, local shell-wall lamina ply- 
gaps, elastic boundary support conditions, and nonuni- 
form loading effects. These results are presented to il- 
lustrate the overall behavior of compression-loaded 
graphite-epoxy shells and the effects of imperfections 
on their response. First, results illustrating the effects 
of orthotropy on the nonlinear response of the 8-ply 
shells are presented. Next, results illustrating the ef- 
fects of initial shell geometric and loading imperfec- 
tions on the buckling response of the 8-ply shells are 
presented. Then, results illustrating the effects of man- 
ufacturing defects in the form of lamina ply-gaps on the 
nonlinear response of the 8-ply quasi-isotropic shell are 
presented. Results illustrating the effects of elastic 
boundary support conditions on the displacement and 
strain response of selected shells are also presented. In 
addition, results that demonstrate the effects of uncer- 
tainties or variations in selected parameters on shell re- 
sponse are presented. These uncertainties include 
uncertainties in the lamina Fiber volume fraction, ap- 
plied loading conditions, and imperfection measure- 
ments. Finally, comparisons between selected 
analytically predicted results and experimentally mea- 
sured results for all four shells are presented. The re- 
sults include predicted and measured load-end- 
shortening response curves, predicted and measured 
load-strain response curves, predicted prebuckling, 
buckling and postbuckling deformation response pat- 
terns, and predicted axial and circumferential stress re- 
sultant patterns. In addition, the results of a local two- 
dimensional stress analysis which identify severe local 
stress gradients in the shell wall of specimen C4, which 
could have caused the specimen to fail before it buck- 
led, are presented. The four composite shells consid- 
ered in this study are interchangeably referred to herein 
as shells or specimens Cl, C2, C3 and C4. 

Effects of Orthotropy and Anisotropy on the Nonlinear 
Response of Geometrically Perfect Shells 

Results from the nonlinear analyses of selected 
geometrically perfect shells are presented in this sec- 
tion. The shells analyzed are based on shells Cl, C2 
and C3 which have orthotropic shell-wall laminates 
[+45/0 2 ] s and [+45/90 2 ] s , and a quasi-isotropic shell- 
wall laminate [+45/0/90] s , respectively. First, the 
load-end-shortening response curves for the shells are 
compared. Next, results illustrating a typical nonlinear 
transient collapse response for the quasi-isotropic shell 
C3 are presented in detail. Finally, results illustrating 
the effects of orthotropy and anisotropy on the nonlin- 
ear response arc presented. Values of the axial load P 
and the end-shortening A presented in this section are 
normalized by the predicted linear bifurcation buckling 


6 

American Institute of Aeronautics and Astronautics 


load for the quasi-isotropic shell C3, - 42,590 

lbs., and the nominal shell-wall thickness, t nom = 0.04 
in., respectively. 

The predicted load-end-shortening responses are 
compared for the three shells in Fig. 9a. The prebuck- 
ling load-end-shortening responses are linear up to the 
general instability point indicated in the figure for each 
shell. The prebuckling slopes of each response curve 
vary according to the individual laminate in-plane axial 
stiffnesses. General instability occurs at normalized ax- 
ial load values of 0.707, 0.868 and 0.979, for shells Cl, 
C2 and C3, respectively. The general instability re- 
sponse is followed by a sudden reduction in the axial 
load supported by the shells which is associated with 
the transient collapse response of the shells. The corre- 
sponding load-time histories of the transient collapse 
response of the shells are shown in Fig. 9b. The initial 
portion of a typical load-time history curve, e.g., the 
load-time history curve for shell Cl, is associated with 
the transition from a state just before buckling occurs 
(time = 0.0) to the state associated with the general in- 
stability point indicated in the figure. This initial por- 
tion of the load-time history curve only occurs in the 
analysis of the geometrically perfect shells considered 
in the present study, and is attributed to the uniformly 
symmetric character of the initial response. The non- 
linear solution algorithm that predicts the structural re- 
sponse of the shell needs a small perturbation in the 
response to be able to change from the axisymmetric 
prebuckling response state to an asymmetric general in- 
stability and transient collapse response state in the 
analysis. Since there are no such geometric or loading 
perturbations in the model for the geometrically perfect 
shells, it is suspected that small numerical variations in 
the transient solution accumulate and provide a small 
perturbation to the equilibrium state of the shell that 
causes the general instability and transient collapse re- 
sponse to occur. Following the initiation of the general 
instability response, a typical load-time history curve 
exhibits a sudden reduction in the axial load supported 
by the shell which is associated with the transient col- 
lapse of the shell. The value of the axial load continues 
to decrease until the axial load attenuates to a steady- 
state value. The kinetic energy in the shell dissipates 
over time and the shell reaches a stable postbuckling 
equilibrium state after approximately 0.035 to 0.04 sec- 
onds have elapsed. Normalized post-collapse post- 
buckling load values equal to 0.403, 0.218 and 0.286 
are obtained for shells Cl, C2 and C3, respectively. 
The results indicate that laminate orthotropy can have a 
significant effect on the load-end-shortening response 
of the shells. In particular, while shell C 1 has the lowest 
general instability load value of all three shells, it also 


has the highest postbuckling load value of all three 
shells. 

The transient deformation responses for selected 
time steps during the transient collapse of quasi-isotro- 
pic shell C3 are presented in Figs. 10a through lOf. The 
deformation responses shown in the figures have been 
scaled for clarity. The shell-wall normal deformations 
at the general instability point shown in Fig. 9 (0.0125 
seconds into the predicted collapse response) change in 
a short period of time from the uniform axisymmetric 
prebucking deformation pattern shown in Fig. 10a to 
the asymmetric transient unstable deformation pattern 
shown in Fig 10b, at which time the shell begins to col- 
lapse. The magnitudes of the shell-wall normal dis- 
placements vary between ±0.5 times the shell-wall 
thickness. The general instability response of the shell 
is caused by the nonlinear coupling of the localized de- 
stabilizing compressive axial and circumferential 
stresses in the shell indicated in Figs. 1 la and 1 lb and 
the normal deformations that occur in the bending 
boundary layer near the ends of the shell. The initial 
buckling deformation pattern indicated in Fig. 10b is 
characterized by localized ellipse-like buckles in the 
bending boundary layer near the ends of the shell. The 
skewing of the deformation response is attributed to the 
presence of the small anisotropic bending-twisting cou- 
pling stiffness terms of the laminate. After 0.013 sec- 
onds have elapsed, the normalized axial load has 
decreased from 0.983 to 0.479, and the deformation re- 
sponse in the shell has evolved to include additional 
rows of small ellipse-like buckles as indicated in Fig. 
10c. The magnitudes of the shell-wall normal displace- 
ments vary from +2 to -4 times the shell-wall thickness. 
After 0.0134 seconds have elapsed, the normalized ax- 
ial load has decreased to 0.465, and the local buckles in 
the deformation pattern have rapidly increased in num- 
ber to include 24 circumferential half-waves and four 
axial half-waves as indicated in Fig. lOd. After 0.0152 
seconds have elapsed, the normalized axial load has de- 
creased further to 0.426, at which time the small buck- 
les in the shell wall begin to coalesce into larger buckles 
with displacement magnitudes that vary from +3 to -7 
times the shell- wall thickness as indicated in Fig. lOe. 
After approximately 0.04 seconds have elapsed, the ki- 
netic energy in the shell has dissipated to a negligible 
level, and the shell has deformed into a stable postbuck- 
ling mode-shape that consists of 16 circumferential 
half-waves and two axial half-waves, as indicated in 
Fig. lOf. 

The transient deformation responses for selected 
time steps during the transient collapse of axially stiff 
orthotropic shell Cl are presented in Figs. 12a through 
12c. The over-all load-time history character of the 
transient collapse response of shell Cl is similar to that 
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of the response exhibited by shell C3, but there are 
some differences in the predicted deformation patterns 
and the collapse initiation events exhibited by shells Cl 
and C3. The deformation pattern of shell C 1 just before 
buckling occurs is characterized by a uniform short 
wave-length axisymmetric response along the entire 
length of the shell as shown in Fig. 12a. These results 
indicate that the bending boundary layer response does 
not attenuate for the laminate of shell Cl, which is in 
contrast to the attenuated bending boundary layer re- 
sponse exhibited by shell C3 as shown in Fig 10a. At 
the general instability point shown in Fig. 9 (after 
0.0167 seconds have elapsed in the transient collapse 
response), a circumferential wave pattern has devel- 
oped into the shell-wall normal deformations for shell 
Cl as indicated in Fig. 12b, which is in contrast to the 
localized response that occurs in the bending boundary 
layer near the ends of shell C3 as shown in Fig. 10b. 
The global collapse of shell Cl initiates in the interior 
of the shell instead of initiating in the bending boundary 
layer as it did for shell C3. In addition, shell Cl does 
not exhibit the pronounced skewing of the deformation 
response that was indicated for shell C3, which sug- 
gests that this shell is less sensitive to the effects of the 
small anisotropic bending-twisting coupling stiffness 
terms of the laminate. After 0.04 seconds have elapsed 
in the transient collapse response, shell Cl exhibits a 
postbuckling deformation pattern that is characterized 
by 16 circumferential half-waves and one axial half- 
wave as indicated in Fig 12c. This pattern is in contrast 
to the pattern of shell C3 which has a postbuckling de- 
formation pattern with 16 circumferential half-waves 
and two axial half-waves. 

The transient deformation responses for selected 
time steps during the transient collapse of circumferen- 
tially stiff orthotropic shell C2 are presented in Figs. 
13a through 13c. The over-all character of the transient 
collapse response of shell C2 is similar to the response 
exhibited by shell C3, however, there are some differ- 
ences in the resulting deformation patterns exhibited by 
shell C2. The deformation pattern for shell C2 just be- 
fore buckling occurs is characterized by a uniform 
short-wave-length axisymmetric bending boundary 
layer response that attenuates rapidly into the interior of 
the shell as shown in Fig. 13a. At the general instabil- 
ity point shown in Fig. 9 (after 0.0217 seconds have 
elapsed in the transient collapse response), the shell ex- 
hibits a deformation pattern that is characterized by one 
axial half-wave and 22 circumferential half-waves as 
indicated in Fig 13b. This pattern is in contrast to the 
short-wave-length responses exhibited by shells C 1 and 
C3 shown in Figs. 12b, and 10b, respectively. In addi- 
tion, the shell-wall deformations exhibit a significant 
amount of skewing that is attributed to the presence of 


the small anisotropic bending-twisting coupling stiff- 
ness terms of the laminate. After 0.037 seconds have 
elapsed in the transient collapse response, shell C2 ex- 
hibits a stable postbuckling deformation pattern that is 
characterized by 14 circumferential half-waves and two 
axial half-wave as indicated in Fig 13c. This pattern is 
in contrast to the pattern for shell C3 which exhibits a 
postbuckling deformation pattern with 16 circumferen- 
tial half-waves and two axial half-waves. 

The previous results indicate that laminate orthot- 
ropy can have an effect on the nonlinear response of 
compression-loaded composite shells. In particular, the 
results indicate that the bending boundary layer attenu- 
ation response can be significantly affected by the 
orthotropy of the laminate. The effects of laminate 
orthotropy on the displacement profiles just before 
buckling occurs are indicated in Fig. 14 for the three 
shells. The dashed, dash-dot, and solid lines represent 
displacement profiles for shells Cl, C2 and C3, respec- 
tively. These results indicate that the axially stiff ortho- 
tropic shell Cl exhibits large-amplitude short-wave- 
length bending deformations just before buckling oc- 
curs. In addition, the bending response decays only 
slightly over the length of this shell. In contrast, shells 
C2 and C3 exhibit short-wave-length bending boundary 
layer responses that attenuate into the interior of the 
shell. These results also indicate that the laminate 
orthotropy can have a significant effect on the magni- 
tude of the maximum displacements in the bending 
boundary layer near the ends of the shell. The displace- 
ments with the largest magnitude are exhibited by shell 
Cl which has a maximum displacement magnitude of 
approximately 0.34 wft nom followed by shells C2 and 
C3 which have maximum displacement magnitudes 
equal to 0.29 and 0.24 wlt nom , respectively. Additional 
results are presented in Fig. 15 that indicate the effects 
of shell length on the bending boundary layer attenua- 
tion characteristics of shell Cl. The solid, dashed, and 
dash-dot lines in the figure represent displacement pro- 
files for shell length-to-radius ratios UR equal to 2, 3 
and 4, respectively. These results indicate that the 
bending boundary layer response for the shell Cl lami- 
nate attenuates for the longer shells, but not for the 
shorter shell. In addition, the results indicate that the 
magnitude of the large-amplitude displacements near 
the ends of the shell decreases as the length of the shell 
increases. These results suggest that the non-attenuat- 
ing bending boundary layers in the shorter shell interact 
to cause displacements with larger magnitudes in this 
shell. These results also suggest that the effects of the 
attenuation length of the bending boundary layer should 
be considered in studies of the nonlinear stability re- 
sponse and failure characteristics of orthotropic shells 
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and in the design of these shells. 

Effects of Initial Imperfections on the Buckling 
Response 

Results from a parametric study that illustrate the 
effects of initial imperfections on the nonlinear re- 
sponse of compression-loaded composite shells are pre- 
sented in this section. The imperfections studied 
include traditional initial geometric shell-wall mid-sur- 
face imperfections and other nontraditional shell-wall 
imperfections and nonuniform load distribution effects. 
The nontraditional shell-wall imperfections include the 
effects of measured shell-wall thickness variations and 
corresponding thickness adjusted lamina properties. 
The nonuniform load distribution conditions include 
the effects of initial shell-end imperfections and mea- 
sured applied load variations. Imperfection-amplitude 
scaling factors or coefficients were introduced into the 
models to allow for independent variations in the im- 
perfection amplitudes during the study. Imperfection- 
amplitude scaling factors for the shell-wall mid-surface 
imperfection, shell-wall thickness imperfection, and 
load distribution imperfection are denoted by W, T and 
D, respectively. Imperfection-amplitude scaling factor 
values equal to 0.0, 0.5 and 1.0 were considered in the 
present study. An imperfection-amplitude scaling fac- 
tor equal to 1 .0 indicates that the actual amplitude oFtKe 
measured imperfection was used in the analysis, and a 
value equal to 0.0 indicates that the imperfection was 
not included in the analysis. An imperfection-ampli- 
tude scaling factor equal to 0.5 indicates that half of the 
measured imperfection amplitude was included in the 
analysis. First, typical results from the parametric study 
illustrating the effects of measured imperfections on the 
buckling load of quasi- isotropic shell C3 are presented. 
Then, typical results illustrating the effects of measured 
imperfections on the nonlinear collapse response of 
shell C3 are presented. 

Results illustrating the effects of measured imper- 
fections on the predicted buckling load of quasi-isotro- 
pic shell C3 are shown in Fig. 16 where the buckling 
load Per is normalized by the linear bifurcation buck- 
ling load P bi f. The open square symbols in the figure 
represent the buckling loads of the shell for different 
values of the parameters considered. The solid and 
dashed lines represent thickness imperfection-ampli- 
tude scaling factors of T equal to 0.0 and 1.0, respec- 
tively. The results indicate that the measured 
imperfections can have a significant effect on the buck- 
ling load of the shell with the normalized buckling 
loads ranging from 0.947 to 0.842. These results indi- 
cate that the magnitude of the buckling load is most sen- 
sitive to the traditional shell-wall mid-surface 
geometric imperfection which can account for as much 


as a 7% reduction in the buckling load of the shell. The 
results also indicate that the shell-wall thickness varia- 
tion can account for a 3-4% reduction in the bucking 
load. In addition, the results indicate that the nonuni- 
form load distribution imperfection can cause a reduc- 
tion in the buckling load of approximately 3-4% for 
small values of the geometric shell-wall mid-surface 
imperfection amplitude. However, as the amplitude of 
the geometric shell-wall mid-surface imperfection in- 
creases, the nonuniform load distribution affects the 
amount of the decrease in the buckling load differently. 
This result indicates that there exists a complex nonlin- 
ear coupling behavior between the shell-wall mid-sur- 
face imperfection and the nonuniform load distribution 
for this shell. In particular, the character of the nonuni- 
form loading is such that, when it couples with the geo- 
metric shell-wall mid-surface imperfections, the shell 
exhibits a redistribution of the internal loads to regions 
of the shell that have higher geometric stiffness values 
which, in turn, retards to onset of buckling in the shell. 
Results for T = 0.5 exhibit similar trends, and the buck- 
ling load values are 1 to 2% lower than those for T = 
0.0. 

Similar results for orthotropic shells C 1 and C2 in- 
dicate that these shells are apparently less sensitive to 
the imperfections considered in the parametric study 
since they exhibit reductions in buckling load values of 
4.1% and 4.7%, respectively, which is in contrast to the 
1 1% reduction exhibited by shell C3. These results also 
indicate that, while the buckling loads are most sensi- 
tive to variations in the traditional shell-wall geometric 
imperfection amplitude, this geometric imperfection re- 
sults in only a 1.7 and a 4.0% reduction in the buckling 
loads of shells Cl and C2, respectively, which is in con- 
trast to the 7% reduction exhibited by shell C3. In ad- 
dition, the shell-wall thickness imperfection and 
nonuniform load distribution imperfection both appear 
to cause no more than a 1% reduction in the buckling 
loads of these shells, which is in contrast to the 3-4% re- 
duction exhibited by shell C3. 

Typical analytically predicted transient response 
deformations that occur during the collapse of shell C3 
with imperfections included in the analysis are shown 
in Figs. 17a through 17d. The imperfection-amplitude 
scaling factors for these results are W= 1 .0, T= 1 .0, and 
D = 1.0. Just before buckling occurs, the shell wall de- 
formations are characterized by several localized el- 
lipse-like deformation patterns as indicated in Fig. 17a. 
The localization in the deformation pattern is caused by 
the combination of a local geometric shell-wall imper- 
fection that is in the form of a significant variation in the 
shell-wall mid-surface geometry, and the intersection 
of a helical ply-gap and a circumferentially aligned ply- 


9 

American Institute of Aeronautics and Astronautics 


gap in the shell at x/L r = 0.25 and 0 = 210°. The local- 
ized deformations occur in regions with destabilizing 
compressive axial and circumferential stresses as indi- 
cated in Fig. 18. After approximately 0.00143 seconds 
have elapsed in the transient response, a single ellipse- 
like buckle has grown in amplitude and couples with 
the destabilizing stresses in the shell wall to cause the 
general instability and collapse of the shell. After ap- 
proximately 0.00238 seconds have elapsed in the tran- 
sient response, additional local buckles have formed 
around the circumference and along the length of the 
shell as indicated in Fig. 17c. As the buckling process 
continues, the deformation pattern in the shell wall con- 
tinues to evolve, and the ellipse-like buckles in the shell 
begin to coalesce into larger diamond-shaped buckles. 
After approximately 0.02 seconds have elapsed in the 
transient response, the kinetic energy in the shell has 
dissipated to a negligible level, and the shell has de- 
formed into a stable postbuckling mode-shape as indi- 
cated in Fig. 17d. These results indicate that the 
collapse response of the imperfect shell is initiated by a 
localized response which leads to the over-all collapse 
of the shell. These results are different from the col- 
lapse response results exhibited by the corresponding 
geometrically perfect shell which are shown in Figs. 
lOa-lOf. The transient collapse response occurs after 
0.00143 seconds have elapsed in the transient analysis 
in the shell with imperfections, and this collapse re- 
sponse occurs much earlier in the load-response history 
of the shell than for the corresponding geometrically 
perfect shell collapse response. This earlier collapse 
time is attributed to the localized prebuckling shell-wall 
deformations and destabilizing in-plane compressive 
stresses in the shell, which result in a rapid transition 
from the stable prebuckling state to the unstable tran- 
sient collapse response. Results for orthotropic shells 
Cl and C2 indicate similar response characteristics as- 
sociated with the transient collapse response of these 
shells. 

Effects of Lamina Plv-Gap Fabrication Defects 

Selected results illustrating the effects of lamina 
ply-gaps on the response of quasi-isotropic shell C3 are 
presented in this section. First, results illustrating the 
effects of lamina ply-gap orientation, ply-gap width, 
and ply-gap depth on the buckling load of this shell are 
presented. Ply-gap orientations of -45°, +45°, 0° and 
90°; ply-gap widths equal to 0. 1 in. and 0.2 in.; and ply- 
gap depths equal to 0.005 in. and 0.0025 in. are consid- 
ered. Then, typical results illustrating the effects of a 
helical ply-gap on the transient collapse response of the 
shell are presented. 

Results illustrating the effects of ply-gap orienta- 


tion, width, and depth on the buckling load of shell C3 
are shown in Fig. 19 where the buckling load P cr is 
normalized by the linear bifurcation buckling load P bi j. 
The open square symbols represent calculated normal- 
ized buckling loads. The solid and dashed lines repre- 
sent result trends for ply-gap depths of 0.005 in. and 
0.0025 in., respectively. The results indicate that the 
buckling load of the quasi-isotropic shell can be affect- 
ed by certain types of lamina ply-gaps. In general, as 
the ply-gap width and ply-gap depth increase, the buck- 
ling load decreases. The results also indicate that the 
buckling load is sensitive to the orientation of the ply- 

gap. These results indicate that shells with 90° ply-gaps 
exhibit the most significant reductions in buckling load 

followed by shells with -45°, +45° and 0° ply-gaps. 
The results also indicate that shells with -45° ply-gaps 
exhibit lower buckling loads than shells with +45° ply- 
gaps. These results suggest that there is a nonlinear 
coupling response that is affected by the laminate 
anisotropy and the local structural response associated 
with the lamina ply-gap. In addition, the results indi- 
cate that the buckling load of the shell is not sensitive to 
the effects of a 0° ply-gap. The benign effect of the 0° 
ply-gap suggests that the local bending deformations 
and destabilizing in-plane stresses near the ply-gap for 
this shell do not significantly affect the buckling of the 
shell. 

Transient deformation response results for select- 
ed time steps during the transient collapse of quasi-iso- 
tropic shell C3 with a -45° lamina ply-gap are presented 
in Fig. 20a through 20d. Just before buckling occurs, 
the shell- wall deformations are characterized by a com- 
bination of a large local inward deformation pattern 
aligned with the ply-gap and a small ellipse-shaped de- 
formation pattern located at the intersection of the 
bending boundary layer deformations near the end of 
the shell and the ply-gap as shown in Fig. 20a. The de- 
stabilizing in-plane axial and circumferential stresses 
shown in Fig. 21a and 21b, respectively, couple with 
the normal shell-wall deformations to cause the general 
instability and collapse of the shell to occur. After ap- 
proximately 0.0012 seconds have elapsed in the tran- 
sient response, additional ellipse-like or diamond- 
shaped buckle patterns have formed in the shell in the 
vicinity of the ply-gap as indicated in Fig. 20b. After 
approximately 0.0029 seconds have elapsed in the tran- 
sient response, additional local buckle patterns have 
formed around the circumference and along the length 
of the shell as indicated in Fig. 20c. As the buckling 
process continues, the deformation pattern in the shell- 
wall continues to evolve, and the ellipsc-Iike buckle 
patterns in the shell begin to coalesce into larger dia- 


10 

American Institute of Aeronautics and Astronautics 


mond-shaped buckle patterns. After approximately 
0.02 seconds have elapsed in the transient response, the 
kinetic energy in the shell has dissipated to a negligible 
level, and the shell has deformed into a stable postbuck- 
ling mode-shape as indicated in Fig. 20d. Similar re- 
sults for orthotropic shells Cl and C2 indicate that the 
nonlinear response and buckling loads of these ortho- 
tropic shells are only slightly affected by the lamina 
ply-gap features. For these orthotropic shells, the ply- 
gaps cause a maximum reduction in the predicted buck- 
ling loads of approximately 1-2%. This benign effect of 
the ply-gaps on the buckling loads of these orthotropic 
shells suggests that the local bending deformations and 
destabilizing in-plane stresses associated with the ply- 
gaps do not significantly affect the buckling loads of 
these shells. 

Effects of Elastic Boundary Conditions 

The potting material applied to the ends of the 
shell specimens is intended to prevent premature dam- 
age to the ends of the shell during testing. In practice, 
it is virtually impossible to achieve perfect rigidly 
clamped supports in a test. It has been shown by Fuchs 
et al. that the displacement and strain responses in cy- 
lindrical shells subjected to an applied bending load can 
be sensitive to the flexibility of the specimen support 
conditions. To determine if the flexibility of the end 
supports may affect the response of compression-load- 
ed shells, an analytical study was conducted to assess 
the effects of axial and radial potting-support stiffness- 
es on the response of the shells considered in the present 
study. Preliminary results indicate that the axial pot- 
ting-support stiffness primarily affects the effective ax- 
ial stiffness of the shell and has no noticeable effect on 
the displacement and strain response of the shell. How- 
ever, the results also indicate that the radial potting-sup- 
port stiffnesses can have a significant effect on the 
displacement and strain response of the shells, and 
some of these results are presented herein. 

Selected results from a parametric study are pre- 
sented to illustrate the effects of elastic radial con- 
straints on the displacement and strain response of 
selected geometrically perfect and imperfect shells. 
Elastic radial stiffnesses K R equal to 1.0E3, 1.0E5and 
1.0E7 lbf/in. were considered in the study. A value of 
the elastic radial stiffness K R - 1.0E5 lbf/in. corre- 
sponds to the effective radial potting-support stiffness 
calculated from the potting model illustrated in Fig. 8, 
and a value of K R = 1.0E7 lbf/in. corresponds to a 
clamped condition. Results illustrating the effects of 
elastic radial potting-support conditions on the preb- 
uckling displacements and strain responses of geomet- 
rically perfect and imperfect orthotropic shell C2 are 


presented subsequently. 

Results illustrating the effects of elastic radial pot- 
ting-support conditions on the typical prebuckling nor- 
mal displacements of the geometrically perfect and 
imperfect orthotropic shell C2 are presented in Figs. 
22a and 22b, respectively. The normal displacement w 
and the axial coordinate x are normalized with respect 
to the nominal shell-wall thickness t nom and the shell 
test-section length L r , respectively. The solid, dashed, 
and dot-dashed lines in the figures represent results for 
shells with elastic radial potting-support stiffnesses K R 
equal to 1.0E7, 1.0E5 and 1.0E3 lbf/in., respectively. 
The results indicate that the variation in the radial 
boundary stiffness can have a measurable effect on the 
magnitude of the displacements for the geometrically 
perfect shell as shown in Fig. 22a. These results indi- 
cate that radial boundary stiffnesses equal to 1.0E5 and 
1.0E7 lbf/in. restrain the radial displacements at the 
boundaries of the shell. However, when the elastic ra- 
dial boundary stiffness is reduced to 1.0E3 lbf/in., the 
boundaries of the shell have an outward radial displace- 
ment of approximately 0.06 w/t nom and a reduction in 
the maximum displacements in the bending boundary 
layer of approximately 0.02 w!t nom . In addition, the 
shell exhibits a slight axial shift in the bending bound- 
ary layer response near the ends of the shell. These re- 
sults indicate that the variation in radial boundary 
stiffness can have a significant effect on the displace- 
ment response of a shell with imperfections as shown in 
Fig. 22b. The results indicate that there can be a 20 to 
25% variation in the local displacement response at sev- 
eral locations along the length of the shell for the differ- 
ent boundary conditions. Other related results indicate 
that these variations in the radial boundary stiffness and 
variations in the prebuckling displacement response 
can cause significant changes in the character of the 
collapse and post-collapse deformation responses of the 
shells. 

Results illustrating the effects of elastic radial 
support conditions on the prebuckling axial strain re- 
sponse of the geometrically perfect and imperfect 
orthotropic shell Cl are presented in Figs. 23a and 23b, 
respectively. The axial strains e^. and the axial coordi- 
nate x are normalized by the axial strain at buckling z cr 
and the shell test-section length Lj , respectively. 

These results indicate that variations in the elastic radial 
support conditions can have an influence on the strain 
response in the geometrically perfect and imperfect 
shells. The results indicate that the strain gradients near 
the ends of the shell shift in the axial direction, which 
cause a local variation in the strain values of as much as 
1 5% and a reduction in the maximum strain of approx- 
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imatcly 5%, In addition, the results indicate that varia- 
tions in the elastic radial support condition can have a 
similar effect on the strain response of an imperfect 
shell as shown in Fig. 23b. These results indicate that 
elastic radial support stiffnesses of 1 .0E3 and 1 .0E5 lbf/ 
in. can cause a local variation in strain values of as 
much as 1 7% near the boundary of the shell. These re- 
sults also suggest that including the effects of elastic ra- 
dial support conditions in the analysis may be important 
when interpreting experimentally measured strain gage 
data for strain gages located near the ends of a speci- 
men. 

Effects of Uncertainties or Variations in Selected Speci- 

Results from a numerical study of the effects of 
uncertainties or variations of selected specimen param- 
eters on the response of the imperfect quasi-isotropic 
shell C3 are presented in this section. The parametric 
uncertainties considered include uncertainties in the 
imperfection measurements, the lamina fiber volume 
fraction, and the applied load distribution. The imper- 
fection measurement uncertainties are attributed to the 
accuracy tolerances of the coordinate measurement de- 
vice used to measure the initial geometry of the shell. 
The shell-wall imperfection measurement values and 
the shell-end-shape imperfection measurement values 
used in the present study are accurate to within ±0.0006 
inches, or ±0.75E-4# and ±1.5E-4L. The shell- wall 
thickness values are accurate to within ±0.0012 inches, 
or ±0.03/ nom . The lamina fiber volume fraction is 
specified by the manufacturer of the material to be 
equal to 0.65 ±0.03 for the 0.005-in. -thick graphite-ep- 
oxy preimpregnated tape material used to fabricate the 
test specimens. Applied load distribution uncertainties 
are measured indirectly by monitoring the measured 
and predicted axial strains at selected points near the 
top and bottom loading surfaces of the shell. A correc- 
tion to the applied displacements can be determined 
from the differences in the measured and predicted 
strains as follows. A user-written program external to 
the STAGS code was used to analyze the differences in 
the measured and predicted strains for a specified ap- 
plied load value. This program used an iterative predic- 
tor-corrector method to determine a correction to the 
applied shell-end displacements. A new finite-element 
analysis was conducted with this displacement correc- 
tion included in the model. This process was repeated 
iteratively until the difference in the measured and pre- 
dicted strains reached a predetermined tolerance. A 
typical predicted shell-end displacement correction for 
shell C3 is shown in Fig. 24. The solid and dashed lines 
represent the predicted top and bottom corrections to 


the displacements applied to the shell, respectively. 
The results indicate that significant corrections to the 
applied end displacements are predicted for this shell, 
and the magnitude of the corrections is on the order of 
the initial shell-end imperfections shown in Fig. 5. The 
majority of this correction is attributed to differences in 
the assumed orientation of the upper loading platen 
with respect to the shell specimen and its actual orien- 
tation prior to applying the load to the specimen. 

Results illustrating the effects of uncertainties in 
se veral specimen parameters on the buckling load of the 
imperfect quasi-isotropic shell C3 are shown in Fig. 25. 
These parametric uncertainties include variations in the 
lamina fiber volume fraction, the measured shell-wall 
thickness imperfection, and the applied load distribu- 
tion. The open square symbols in the figure represent 
calculated normalized buckling loads, where the buck- 
ling load P cr is normalized by the linear bifurcation 
buckling load P^j. The solid, dashed, and dot-dashed 
lines in the figure represent trends in the results for 
shells with thickness imperfection measurement varia- 
tions of -3%, 0%, and +3%, respectively. These results 
indicate that the normalized buckling loads range from 
0.83 to 1.0, and that variations in these parameters can 
have a significant effect on the buckling load of the 
shell. In particular, the results indicate that a variation 
in the shell wall thickness of ±3% can result in a 6% 
variation in the buckling load of the shell. Similarly, re- 
sults indicate that a ±3% variation in the fiber volume 
fraction can result in an 8% variation in the buckling 
load. In addition, the buckling load is decreased by up 
to 2% when the magnitude of the applied load distribu- 
tion correction amplitude scaling factor is increased 
from 0.0 to 1.0. Results for shells Cl and C2 indicate 
similar tends. However, these results suggest that the 
orthotropic shells are less sensitive to predicted applied 
load distribution uncertainties than is the quasi-isotro- 
pic shell. 

Predicted and Measured Response Com parisons 

Selected results from nonlinear analyses of the 
orthotropic and quasi-isotropic shells are compared to 
the experimentally measured results in this section. 
The nonlinear analysis results are for shell models that 
included the effects of the measured initial shell-wall 
geometric and thickness imperfections, thickness-ad- 
justed material property variations, shell-wall lamina 
ply-gaps, measured loading variations, elastic radial 
support conditions, and selected specimen parametric 
uncertainties. A value of the elastic radial support stiff- 
ness equal to 1.0E5 Ibf/in. is used in the models to sim- 
ulate the boundary conditions provided by the potted- 
end supports of the shell specimens. Upper and lower 
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response bounds were determined based upon the re- 
sults of a traditional combinatorial analysis of the ef- 
fects of selected specimen parametric uncertainties. 
Uncertainties in the measured thickness imperfection, 
lamina fiber volume fraction, and applied load distribu- 
tion were considered in the analysis. Predicted and 
measured load-end-shortening response curves, select- 
ed load-strain responses, and observed and predicted 
post-collapse normal displacements are presented in 
this section. 

Analytically predicted and experimentally mea- 
sured load-end-shortening response curves for shells 
Cl, C2 and C3 are presented in Fig. 26, Values of the 
axial load P and the end-shortening A are normalized 
by the linear bifurcation buckling load of quasi-isotro- 
pic shell C3, P = 42,590 lbs., and the nominal shell- 
wall thickness t nom - 0.04 in., respectively. The solid 
and dashed lines in the figure represent experimentally 
measured and analytically predicted results, respective- 
ly. Each shell has two predicted response curves repre- 
senting analytically predicted upper and lower bounds 
to the response based on specimen parametric uncer- 
tainties. The measured results indicate that shells Cl, 
C2 and C3 exhibit the general instability points shown 
in the Figure at normalized load values of 0.652, 0.749 
and 0.803, respectively, and are 7.8, 13.7 and 17.6% 
lower than the predicted loads for the geometrically 
perfect nominal shells shown in Fig. 9a, respectively. 
These results indicate that, for the most part, the mea- 
sured load-end-shortening response curves for the 
shells fall on or within the analytically predicted re- 
sponse bounds. These results indicate that the mea- 
sured response curves tend to correlate well with the 
analytically predicted lower bounds for each shell. An- 
alytically predicted and experimentally measured strain 
response curves for selected points on shell C3 are pre- 
sented in Figs. 27a and 27b. Values of the axial load P 
and the axial strain t x are normalized by the linear bi- 
furcation buckling load of the quasi-isotropic shell C3, 
Pty* “ 42,590 lbs., and the corresponding axial strain 
at buckling e CT = 0.00257, respectively. The strain re- 
sults in both figures represent back-to-back surface 
strain measurements within the bending boundary layer 
near the ends of the shell. Strain gages 79 and 80 are 
outer and inner surface strain gages located at x/L j- 

0.071 and 0 = 0°, and strain gages 83 and 84 are outer 
and inner surface strain gages located at x/L r = 0.071 

and 9 = 180° The bending boundary layer behavior is 
evident in the divergent character of the back-to-back 
predicted and measured strain results. These results in- 
dicate that the measured prebuckling strain responses, 
for the most part, fall on or within the predicted re- 


sponse bounds. However, there are some discrepancies 
in the predicted and measured post-collapse strain re- 
sponses that are discussed subsequently. The results 
also indicate that the predicted post-collapse strain re- 
sponse can be very sensitive to specimen parametric un- 
certainties as shown in Fig. 27b. The measured and 
predicted prebuckling strain results for orthotropic 
shells Cl and C2 agree well overall, but there are some 
discrepancies in the post-collapse predictions. 

Predicted initial post-collapse normal displace- 
ment contours and the corresponding observed moird 
fringe patterns for specimen Cl are shown in Fig. 28. 
The dashed contour lines in the predicted displacement 
contour plots represent inward displacements and the 
solid lines represent outward displacements. The den- 
sity of the contour lines indicates the severity of the dis- 
placement gradients in the specimen. These results 
indicate that the specimen collapses into a general-in- 
stability diamond-shaped buckling pattern with 16 half- 
waves around the circumference and one half-wave 
along the length, as predicted by the transient analysis. 
The analytical results predict that the mode-shape pat- 
tern is 15-20° out-of-phase with the observed mode- 
shape pattern. Predicted initial post-collapse normal 
displacement contours and the corresponding observed 
muir£ fringe patterns for specimen C2 are shown in Fig. 
29. These results indicate that the specimen collapses 
into a general instability diamond-shaped buckling pat- 
tern with 14 half-waves around the circumference and 
two half-waves along the length, as predicted by the 
transient analysis. In addition, the predicted mode- 
shape is in-phase with the observed mode-shape. Sim- 
ilar results for shell C3 indicate that the shell collapses 
into a general-instability diamond-shaped pattern with 
16 circumferential half-waves and two axial half- 
waves. However, the analytical results predicted that 

the mode-shape pattern is approximately 15° out-of- 
phase with the observed moir£ fringe pattern. This dis- 
crepancy explains, in part, the large discrepancies in the 
post-collapse strain response indicated in Fig. 27. 

Analytically predicted and experimentally mea- 
sured load-end-shortening response curves for the 8- 
ply quasi-isotropic shell C3 and the 16-ply quasi-isotro- 
pic shell C4 are shown in Fig. 30. The axial load P is 
normalized by the quantity EA, i.e., the effective axial 
stiffness of the shell, denoted by E, multiplied by the 
nominal shell cross-sectional area, denoted by A , and 
the end-shortening A is normalized by the nominal 
shell length L= 16.0 in. The solid and dashed lines in 
the figure represent experimentally measured and ana- 
lytically predicted results, respectively. Each shell has 
two predicted response curves representing analytically 
predicted upper and lower bounds to the response based 
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on specimen parametric uncertainties. The analytical 
results predict that shell C4 will buckle at an effective 
axial strain value (AIL) of approximately 0,0052- 
0.0054, which is approximately twice the value predict- 
ed for shell C3. The results for shell C3 indicate that for 
the most part, the measured load-end-shortening re- 
sponse curve for shell C3 fall on or within the analyti- 
cally predicted response bound, however. The 
measured results for shell C4 fall outside of the predict- 
ed bounds for effective strain values greater than ap- 
proximately 0.003, and the shell has a limit load that is 
10% lower than the predicted lower bound buckling 
load for the shell. The measured limit load for shell C4 
corresponds to the complete failure of the shell and is 
characterized by significant material damage around 
the circumference of the shell. This shell had no post- 
failure residual strength. Post-test inspection of the 
shell indicated that the overall failure of the shell may 
have been initiated by a material failure near a circum- 
ferentially aligned ply-gap in the shell wall. To inves- 
tigate that character of a possible material failure, a 
two-dimensional generalized plane-strain model of the 
ply-gap feature was developed to predict the local stress 
distribution around the ply-gap. The ply-gap model ge- 
ometry and loading conditions are shown in Fig. 31a. 
The local ply-gap feature is subjected to a combined 
loading condition which is predicted by the global fi- 
nite-element analysis of the entire shell. The loads ap- 
plied to the ply-gap model are representative of the 
local loads near the ply-gap at buckling. The applied 
loads include an axial strain z x equal to 0.005, a bend- 
ing rotation r y equal to 0.0002 radians, and a trans- 
verse shear load Q x equal to 20.0 lbf/in. The predicted 
axial normal stress <7 X , transverse normal stress o z , and 
transverse shear stress c xz from the analysis are shown 
in Fig 31b, c and d, respectively. The results indicate 
that the ply-gap causes significant stress concentrations 
to develop within the laminate. More specifically, the 
0° lamina ply near the ply-gap and the 0° ply that exists 
towards the outer surface of the laminate have elevated 
values of all the stress components. The predicted re- 
sults indicate that stress levels are equal to or exceed the 
typically reported stress allowables for this material 
system, and suggest that fiber compression failures and 
delamination failures could have occurred for these ap- 
plied load values. In particular, several regions of the 
laminate have transverse shear stress values that are on 
the order of two to three times the allowable shear stress 
value for the material, and these results suggest that 
delamination failures near the ply-gap may have oc- 
curred for load values less than the predicted buckling 
load of the shell. These results suggest an explanation 
for some of the discrepancy between the predicted and 


measured results for shell C4. 

Design Considerations 

The nonlinear analysis procedure described herein 
offers a robust and accurate approach for predicting the 
nonlinear response and stability characteristics of com- 
pression-loaded thin-walled composite shell structures. 
This nonlinear analysis procedure can be used to form 
the basis of a modern high-fidelity design and analysis 
approach for composite shell structures that accounts 
for the effects of imperfections associated with the fab- 
rication process used to fabricate the structure. A com- 
parison of results from the traditional shell design 
approach and from the present nonlinear high-fidelity 
shell analysis approach is shown in Fig. 32 for shell C3. 
Values of the axial load P and the end-shortening A 
in the figure are normalized by the predicted linear bi- 
furcation buckling load for quasi-isotropic shell C3, 

P q ™ 51 = 42,590 lbs., and the nominal shell-wall thick- 
ness t nom = 0.04 in., respectively. The traditional ap- 
proach to shell design is to predict the shell buckling 
load using a linear bifurcation buckling analysis with 
the nominal structural dimensions and material proper- 
ties of an idealized geometrically prefect shell. This 
predicted bucking load for the shell is then reduced by 
an empirical “knockdown” factor based on a design cri- 
terion such as the lower-bound design recommenda- 
tions reported in NASA SP-8007 (Ref. 1). Quasi- 
isotropic shell C3 has a nominal radius-to-shell-wall- 
thickness ratio equal to 200, and the design knockdown 
factor from NASA SP-8007 for an equivalent isotropic 
shell is approximately equal to 0.47 as indicated in Fig. 
32. This traditional design approach results in an overly 
conservative design for this shell. The normalized mea- 
sured buckling load and the normalized predicted buck- 
ling load from the lower response bound predicted by 
the nonlinear analysis procedure are 0.803 and 0.832, 
respectively, and these results are approximately 33 and 
36% greater than the results based on the traditional de- 
sign approach. The modern high-fidelity analysis ap- 
proach offers a relatively affordable alternative to 
relying on historical test data for shells that do not rep- 
resent the configuration, material system, or fabrication 
process for a particular composite shell design of inter- 
est. This suggested design analysis procedure should 
be used with a selected number of carefully conducted 
experiments that could be used to verify the design and 
analysis results. This approach could make it possible 
to avoid testing the large number of replicates of a de- 
sign needed to develop empirical design factors. This 
nonlinear analysis procedure could be used as a para- 
metric tool in the early stages of a design development 
program to determine the sensitivity of the response 


14 

American Institute of Aeronautics and Astronautics 


program to determine the sensitivity of the response 
characteristics of a specific design to a number of dif- 
ferent types of imperfections or differences in the ideal- 
ized as-designed shell structure and the actual as- 
manufactured shell structure. Since this nonlinear anal- 
ysis procedure can predict the local stresses and strains 
in a shell for any point in the load-response history of 
the shell, failure of the shell due to local stress or strain 
gradients could also be predicted. 

Concluding Remarks 

The results of an experimental and analytical 
study of the effects of imperfections on the nonlinear 
response and buckling loads of unstiffened thin-walled 
compression-loaded graphite-epoxy cylindrical shells 
with four shell-wall laminates are presented. The shell- 
wall laminates considered in this study include two qua- 
si-isotropic laminate and two different orthotropic lam- 
inates. Shell-radius-to-thickness ratios equal to 100 and 
200 were also considered. The results for the nonlinear 
prebuckling, buckling, transient collapse, and post-col- 
lapse response of geometrically perfect shells and shells 
with measured imperfections are presented. The results 
identify the effects of traditional initial geometric shell- 
wall mid-surface imperfections and the effects of other 
nontraditional imperfections on the nonlinear response 
and buckling loads of the shells. These nontraditional 
imperfections include shell-wall thickness variations, 
material property variations, shell-end geometric im- 
perfections, local shell-wall ply-gaps associated with 
the fabrication process, variations in loads applied to 
the end of the shell, and elastic boundary support con- 
ditions. In addition, upper and lower bounds to the non- 
linear response of the shells are presented which were 
determined from a combinatorial analysis of the effects 
of uncertainties or variations in several shell parame- 
ters. A high-fidelity nonlinear shell analysis procedure 
has been used to predict the nonlinear response of the 
shells, and the analysis procedure accurately accounts 
for the effects of these traditional and nontraditional im- 
perfections and elastic boundary conditions on the non- 
linear response and buckling loads of the shells. The 
analysis results generally correlate well with the exper- 
imental results indicating that it is possible to predict 
accurately the complex nonlinear response and buck- 
ling loads for compression-loaded composite shell 
structures when the traditional and nontraditional im- 
perfections considered in the present study are included 
in the analysis. 

The analytical results indicate that the effects of 
the traditional and nontraditional imperfections consid- 
ered in this study can be important for predicting the 
buckling loads of composite shells since they can sig- 
nificantly affect the nonlinear response and buckling 


loads of the shells. The results indicate that the mea- 
sured imperfections can couple with the in-p!ane com- 
pressive stress resultants in a nonlinear manner to affect 
the shell response. The results show that shell-wall 
laminate orthotropy can also have a significant effect on 
the nonlinear response and buckling loads of composite 
shell structures. However, the results indicate that the 
nonlinear response and buckling loads of the orthotro- 
pic shells studied are less sensitive to the measured im- 
perfections considered in this study than are the 
nonlinear response and buckling loads of the quasi-iso- 
tropic shell studied. In addition, the results indicate that 
elastic boundary support conditions can have a signifi- 
cant effect on the displacement and strain responses of 
the shells; shell length and shell-wall orthotropy can 
have a significant effect on the attenuation of the shell- 
end bending boundary layer; and shell-wall ply-gaps 
associated with the fabrication process can have a sig- 
nificant effect on the shell buckling loads. Moreover, 
the results indicate that shell-wall ply-gaps can cause 
significant stress gradients to develop in the shell that, 
in turn, could cause premature failure of the shell, 

The results indicate that, for the most part, the 
measured response of the shells falls on or within the 
predicted upper and lower bounds to the response that 
are associated with the uncertainties or variations in the 
shell parameters considered in the study. These results 
indicate that the nonlinear analysis procedure used in 
this study can be used to determine accurate, high-fidel- 
ity design knockdown factors that can be used for pre- 
dicting composite shell buckling and collapse loads in 
the design process. The traditional and nontraditional 
imperfections considered in this study could be used to 
formulate the basis for a generalized imperfection sig- 
nature of a composite shell that includes the effects 
variations or uncertainties in the shell-geometry, fabri- 
cation-process, load-distribution and boundary stiffness 
parameters. The high-fidelity nonlinear analysis proce- 
dure used in this study can be used to form the basis for 
a shell analysis and design approach that includes this 
generalized imperfection signature and addresses some 
of the critical shell-buckling design criteria and design 
considerations for composite shell structures without 
resorting to the traditional empirical shell design ap- 
proach that can lead to overly conservative designs. 
Since the nonlinear analysis procedure can be used to 
predict local shell-wall stresses and strains at any point 
in the shell load-response history, the analysis proce- 
dure can also be used to form a robust failure analysis 
for composite shell structures with nonlinear response 
characteristics. 
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Fig. 3 Typical measured wall thickness variation for a shell specimen C3. 
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Fig. 5 Topical measured shell-end or loading-surface imperfections for a shell specimen C3. 
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Fig. 6 Typical experimentally measured loading platen rotations. 
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Fig. 7 Idealized ply-gap model detail. 
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Fig. 8 Topical finite-element model and geometry and boundary conditions of the potting-shell detail and typ- 
ical predicted axial strain distribution at the mid-surface of the shell wall. 
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Fig. 9 Numerically predicted nonlinear response of geometrically perfect compression-loaded cylindrical 
shells. 



a) Prebuckling deformations 
Time = 0.0 seconds 
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b) Initial collapse deformations 
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c) Transient collapse deformations 
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d) Transient collapse deformations 
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Fig. 10 Numerically predicted collapse response for a geometrically perfect compression-loaded quasi-isotro- 
pic shell C3. 
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Fig. 11 Numerically predicted axial and circumferential stress resultants, N x and N r respectively, just before 
buckling for a geometrically perfect compression-loaded quasi-isotropic shell C3. 
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Fig. 12 Numerically predicted collapse response for a geometrically perfect compression-loaded orthotropic 
shell Cl. 



a) Prebuckling deformations 
Time = 0.0 seconds 
P/P, 7" - 0.868 


b) Initial collapse deformations 
Time = 0.02 1 7 seconds 

PIP*™' = 0.866 


c) Post-collapse 
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Fig. 13 Numerically predicted collapse response for a geometrically perfect compression-loaded orthotropic 
shell C2. 
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Fig. 14 Effects of orthotropy on displacements just before buckling of geometrically perfect compression- 
loaded composite cylindrical shells. 



Fig. 15 Effects of shell length on bending boundary layer attenuation for a geometrically perfect compression- 
loaded orthotropic shell Cl. 
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Fig. 16 Effects of measured imperfections on the normalized buckling load of a compression-loaded quasi-iso- 
tropic cylindrical shell C3. Imperfection amplitude scale factors for the measured geometric shell-wall mid- 
surface imperfection, shell-wall thickness imperfection, and nonuniform loading are denoted by W, T> and Z), 
respectively. 
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Fig. 17 Numerically predicted collapse response for an imperfect compression-loaded quasi-isotropic shell C3. 
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a) Axial stress resultant contours, N x , M/in, 
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b) Circumferential stress resultant contours, N y , Ibf/in. 


Fig. 18 Numerically predicted axial and circumferential stress resultants, N x and N y respectively, just before 
buckling for an imperfect compression-loaded quasi-isotropic shell C3. 
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Fig. 19 Effects of ply-gap angle, gap width gw> and gap depth gd , on the normalized buckling load of a geomet- 
rically perfect compression-loaded quasi-isotropic cylindrical shell C3. 
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Fig. 20 Effects of a 45° lamina ply-gap on the transient collapse response of a compression-loaded quasi-iso- 
tropic cylindrical shell C3. 
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a) Axial stress resultant contours, N x , lbf/in. 
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b) Circumferential stress resultant contours, N y , lbf/in. 


Fig. 21 Effects of a 45° lamina ply-gap on the axial and circumferential stress resultants incipient to buckling 
for a compression-loaded quasi-isotropic cylindrical shell C3. 
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Fig. 22 Effects of elastic radial support conditions on the buckling displacement profile of compression-loaded 
orthotropic shell C2. 




Fig. 23 Effects of elastic radial support conditions on the axial strain distribution of a compression-loaded 
orthotropic shell Cl. 
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Fig. 24 Typical numerically predicted applied displacement corrections for the top (8^ ) and bottom (8^J ) 
specimen ends that are inferred from experimentally measured axial strains for quasi-isotropic shell C3. 



Fig. 25 Effects of variations in fiber volume fraction, measured shell-wall thickness, and applied load on the 
normalized buckling load of quasi-isotropic shell C3. 
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Fig. 26 Numencally predicted and experimentally measured Ioad-end-shortening response curves; predicted 
results represent response bounds. 



Predicted 

bounds 


Measured 

result 


P 

y quasi 0.5 


. __ A Strain gage 80 

Strain gage 79 S \ y y 


L V* 

0.75 k 


Predicted 

bounds 


Measured 

result 


vV^Strain gage 84 


Strain gage 83- 


Fig. 27 Numerically predicted and experimentally measured load-strain response curves for quasi-isotropic 
shell C3; predicted results represent response bounds. 
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a) Observed post-collapse moire fringe pattern. 


b) Numerically predicted post-collapse normal 
deformation response. 

Fig. 28 Observed and predicted initial post-collapse normal displacements for specimen Cl. 


220 ° 220 ° 



a) Observed post-collapse moire fringe pattern. b) Numerically predicted post-collapse normal 


deformation response. 

Fig. 29 Observed and predicted initial post-collapse normal displacements for specimen C2. 
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c) Transverse stress contours, <T Z , Ibf/in 2 . 


d) Transverse shear stress contours, a xz , Ibf/in 2 . 


Fig. 31 Typical two-dimensional ply-gap model and numerically predicted through-the-thickness stress con- 
tours near a ply-gap in the sixteen-ply quasi-isotropic shell C4. 
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Fig. 32 Numerically predicted and experimentally measured load-strain response curves for quasi-isotropic 
shell C3. 
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